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The existence of a generalized fluctuation-dissipation theorem observed in simulations and exper- 
iments performed in various glassy materials is related to the concepts of local equilibration and 
heterogeneity in space. Assuming the existence of a dynamic coherence length scale up to which 
the system is locally equilibrated, we extend previous generalizations of the FDT relating static 
to dynamic quantities to the physically relevant domain where asymptotic limits of large times 
and sizes are not reached. The formulation relies on a simple scaling argument and has thus not 
the character of a theorem. Extensive numerical simulations support this proposition. Our results 
quite generally apply to systems with slow dynamics, independently of the space dimensionality, the 
chosen dynamics, or the presence of disorder. 
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How slow the wind 
How slow the sea 
E. Dickinson 

Materials can be far from their equilibrium state for 
two main reasons. First, the relaxation time can be very 
large compared to the experimental time scale and the 
system does not reach equilibrium. This is typically the 
case for a liquid supercooled through its melting transi- 
tion below its glass transition. The glass state is thus a 
non-stationary non-equilibrium state, where the glass is 
said to age. 'Glass phases' are extremely commonly en- 
countered in condensed matter pj. Second, when some 
external drive mechanism is applied to the system, it also 
leaves its equilibrium state. Typical situations can be 
driven interfaces, or soft materials perturbed by a me- 
chanical flow. In this paper, we consider these two types 
of situations in the specific case of systems with slow 
dynamics. We discuss theoretically and numerically the 
existence and formulation of a generalization of the equi- 
librium fluctuation-dissipation theorem and its link to 
equilibrium properties. 

Important results for the understanding of the non- 
equilibrium dynamics of glassy materials have been ob- 
tained in the last decade from the asymptotic solution 
of the dynamics of various infinite-range glassy models 
in the two types of situations described above Q. It 
is established that much information on the dynamics 
is gained by studying two-time correlation and response 
functions, the relationship between them, and the link 
between dynamic and static properties. In the case of an 
aging dynamics, for instance, two-time functions explic- 
itly retain a dependence on their two time arguments, 
and conjugated correlation and response functions do 
not satisfy the relation they follow at equilibrium, the 
fluctuation-dissipation theorem (FDT). In the case of a 
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driven dynamics, time translation invariance might be 
preserved even in the glassy phase, but the equilibrium 
FDT is not satisfied. 

These 'mean-field' results have been used as a guide to 
interpret numerical and experimental studies of realistic 
glassy materials 0- However, since no exact solution is 
known in that case, a naive mean-field interpretation of 
experimental results can be misleading. It is our aim to 
discuss the interpretation of experimental and numerical 
measurements of correlation and response functions in 
terms of static quantities. For this purpose, we extend 
previous analysis to the physically relevant domain where 
large times and sizes are not accessible. Our approach, 
which relics on the concepts of a spatially local equilibra- 
tion and heterogeneity of the system, closely follows the 
line of thought suggested in Refs. 0, H SH H H E3, 
where the need to take into account spatial aspects of 
slow dynamics was already emphasized in this particular 
context. 

The paper is organized as follows. In Section [IJ we 
define the static and dynamics quantities of interest and 
briefly review the known relevant results. In Section ITT1 
we use a scaling argument to support our generalization 
of the FDT. Numerical simulations of aging and driven 
dynamics in a spin glass model in finite spatial dimen- 
sions d = 3 and d — 4 are reported in Section ITTT1 Sec- 
tion ]W\ concludes the paper. 



I. GENERALIZED FDT: WHAT DO WE KNOW? 

In this section, we introduce the basic concepts and 
quantities of interest for this paper. The exact results ob- 
tained in infinite-range glass models are first presented. 
We then review the argument given to support a possi- 
ble larger validity of the mean-field results, and conclude 
the section with the experimental and numerical evidence 
available. The goal of this section is also to recall what 
is known, what is only supposed, and most importantly 
to emphasize what is not known. 
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A. Infinite-range glass models 

The solution of the non-equilibrium dynamics of 
infinite-range glass models consists of the asymptotic 
analysis of coupled dynamical equations involving two- 
time correlation functions, C(t\, t 2 ), and their thermody- 
namically conjugated response functions, R{t\,t 2 ) 
Alternatively, integrated responses, or susceptibilities, 
Xylite) = Jj 1 dt'R(tx,t'), can be studied, with the ad- 
vantage that ^2) is easier to measure in a simulation 
or an experiment. Adopting notations appropriate to 
magnetic systems, those functions read for instance 
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C{tx,t 2 ) = jy2*i(tl)*i(t2 
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h=0 



for a system composed of N spins Si(i = 1, • • ■ , N). The 
fields hi are thermodynamically conjugated to the spins. 
Since and (0) are self-averaging, no average of any 
kind is required if the thermodynamic limit, N — > 00, is 
taken. 

It turns out that in infinite-range glass models, a gen- 
eralized form of the FDT is satisfied in the asymptotic 
limit of large times in the aging case, or small driving 
forces for the driven case. The limits respectively read 
t\, t 2 — * co with C(t\, t 2 ) fixed, or e — > with C(t\ — t 2 ) 
fixed, where e generically refers to the amplitude of the 
driving force (see below for examples). This asymptotic 
form of the FDT reads 



R(h,t 2 ) = 



X[C{t lt t2)] dC{hM) 



T 



dt 2 



(3) 



where T is the temperature, and X(C) is called the 
fluctuation-dissipation ratio (FDR). The equilibrium 
FDT is recovered when X(C) = 1. The finding of a 
generalized FDT for non-equilibrated systems suggests 
the possibility to develop a non-equilibrium extension of 
standard statistical mechanics or thermodynamics. In 
particular, the quantity T/X has been shown to play the 
role of a non-equilibrium effective temperature for the 
slow modes of the system [l^, leading to an important 
body of related works . 

Moreover, in infinite-range glass models which stati- 
cally exhibit a full replica symmetry breaking pattern, 
the new dynamical quantity X(C) was found to be di- 
rectly related to the static order parameter involved in 
the replica symmetry breaking, the distribution of over- 
laps P(q), through the very simple relation [Tl| . 



X(C)= [ dg'P(g') 
Jo 



The 'Parisi function' is defined as 



P(q) = lim 

N~>oc 



N ^ 



(4) 



(5) 



where sf denotes the value of spin i in the configuration 
a, and (• • •) represents an average over configurations 
(a, (3) weighted by their Boltzmann probability. 

Equation has a remarkable form, since it relates two 
quantities of completely different origin. The FDR in the 
left hand side of 0) is measured in the off-equilibrium 
dynamics, while the Parisi function in the right hand 
side is computed statically using the equilibrium Gibbs 
measure. Note also that the computation of X(C) does 
not involve replicas. 

Finally, we recall the remarkable result that a gener- 
alized FDT as in Eq. © is asympt otically obeyed in the 
case of the driven dynamics [13, uM > with the same value 
of the FDR, so that Eq. Q remains unchanged in that 
case Ha. 



B. Stochastic stability 

Further developments have taken a more speculative 
form, in the sense that no exact result has been obtained 
for finite dimensional models. However, it was argued 
that Eq. ffl is generally true for finite dimensional glassy 
systems [3J . We now briefly review the formal derivation 
of this result to emphasize the main hypothesis which are 
made. Following Refs. [TT|| . the generalized susceptibili- 
ties x r are fi rs t introduced: 



X' 



ii<-<i T 



dhi 



(6) 



h=0 



The 'derivation' of Eq. (@J simply consists of the com- 
parison between a static and a dynamic computation of 
X r 0- A static average gives 



( X r ) stat = lim lim ( X r (t)) = i r&qP{q){q r I). 

N^oc t— *oo ± Jq 



A dynamical average leads instead to 



(7) 



(X r )dy„ = lim lim ( X r (t)) = i 

t— >co N^oo 1 



(I 



dC 



Equation directly follows by requiring the equality 
(X r )stat = (x r )dyn for all r. 

There are, however, several assumptions in this deriva- 
tion. First, one assumes that a non-trivial form of the 
FDT given by Eq. J2J is valid in the large time limit. Sim- 
ilarly, the static computation is performed in the thermo- 
dynamic limit, N — > 00. A third assumption is that static 
and dynamic calculations give the same result. This can 
be justified by a standard nucleation argument 0], al- 
though exceptions are known [llj . The last and strongest 
hypothesis has been called 'stochastic stability' 3J. The 
name originates from the fact that the susceptibility © 
probes the behaviour of the system under the random 
perturbation 



pcrturbatii 



= 5 



E ■ 

il<---<ir 



(9) 
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in the limit S — > 0. Here, the field h^...^ is a random 
Gaussian variable of mean and variance r\/(2N r ~ 1 ). 
Stochastic stability implies thus that the limits 5 — > 
[or equivalently h — > in Eq. ©], TV — > oo, and t — > oo 
can all be inverted. This last assumption is completely 
uncontrolled, so that Eq. J3J is made physically plausible 
by this calculation, but has no 'theorem' character. A 
measurable test of this last hypothesis in a realistic sys- 
tem is to check the relation Q it implies, thus making 
the argument a bit circular. Interesting examples and 
counterexamples were discussed in detail in Ref. • 



C. Numerical and experimental results 

Since the generalization of the FDT through Eqs. © 
and 10} cannot be rigorously proved, numerical simula- 
tions and experiments are necessary in order to make 
some progress. The advantage of simulations is that 
both equations can be separately checked, while exper- 
iments can obviously not compare dynamical to static 
data, since there are no experimental way to directly 
measure the Parisi function. Its indirect determination 
through Eq. (QJ would thus be a major result, as pointed 
out in Ref. [3|. It is the principal aim of this paper to 
discuss the issue of what is precisely determined from a 
dynamical measurement performed in a physically acces- 
sible time window. 

Dynamic results have been numerically obtained in 
many different glassy models, as reviewed in Ref. [l3| . 
More recently, experiments have also been able to in- 
dependently probe two-time susceptibilities and correla- 
tions in a number of glassy materials 0, 0, 0, Il9| . 
The standard way to present data and test Eq. is to 
build an 'FD plot' of the susceptibility x(ti> *2) as a func- 
tion of the correlation function C(ti,t2), parameterized 
by the time difference t\ — t%, conventionally taken as 
positive H3- If Eq. © holds, one gets indeed 



X(*l>*2) 



dqX(q), 



(10) 



C(ti,t 2 ) 



implying that the x(C) relation is independent of time, 
and has a slope related to the FDR, 



dx 
dC 



X{C) 
T ' 



(11) 



We have implicitly assumed that C{t\,ti) — 1 to simplify 
the notations. It is important to notice that neither simu- 
lations nor experiments have ever reported an asymptotic 
FD plot, i.e., a time independent relation between x an d 
C . Quite strikingly then, physically accessible time scales 
are such that Eq. © is not valid. As a consequence, it is 
necessary to theoretically investigate preasymptotic be- 
haviours where times are large, but finite. Equivalently 
for the driven dynamics, the case of small but finite driv- 
ing forces has to be considered. 



Far less works have investigated the validity of Eq. 0} . 
Numerically, it is possible to compute the Parisi function 
in an equilibrium simulation, obtaining thus the function 
P(q,L), where L is the — necessarily finite — linear size of 
the system. The L-dependence is explicitly kept since it 
is again important to note that no numerical work has yet 
been able to report the asymptotic form of this function, 
P(q,L — > oo), in a glassy system. Moreover, the system 
size will play an important role below. 

Early studies have then empirically compared static 
and dynamic data through the comparison of the curves 
x(C, £2) obtained in the dynamics for times ti 'as large as 
possible', to the static curve S(C, L) measured statically 
in a system of size L 'as large as possible' [2(3, |^|. We 
have defined 



S(C,L) = ± J c &lf dq'P(q',L). 



(12) 



With these definitions, Eq. 0} can be synthetically 
rewritten as 



X (C) = S(C), 



(13) 



where S(C) = lim^oo S(C, L). 

The remarkable coincidence between static and dy- 
namic data initially found in numerical simulations of 
finite dimensional spin glasses in dimensions d — 3 and 
d = 4 was first taken as the sign that both quantities had 
converged to their asymptotic limits |20j . Conclusions on 
the correct description of the low-t emp erature phase in 
spin glasses phase were then drawn |20| . 

These conclusions were however premature, as strik- 
ingly demonstrated in recent experiments which have 
shown that even in the experimental time window the 
FD plots still have a clear time dependence This 
again emphasizes the necessity of a careful theoretical 
study of preasymptotic behaviours in order to get a cor- 
rect interpretation of experimental data. This issue is 
discussed throughout the rest of the paper. 



II. GENERALIZED FDT WITH SPACE 

A. What about space? A scaling argument 

To understand preasymptotic behaviours, a physical 
picture of the slow dynamical processes involved in the 
system is needed, since infinite-range glass models pro- 
vide us with no relevant prediction in that regime. Un- 
fortunately, due to their mean-field nature, they do not 
provide us with a clear description of the physics either. 
One has thus to resort to more phenomenological descrip- 
tions. 

It is now well established, at least in the case of spin 
glasses, that slow dynamics is accompanied by the ex- 
istence of a 'coherence length' scale, £: the larger the 
coherence length the slower the dynamics. Here, an anal- 
ogy with coarsening phenomena in non-disordered ferro- 
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magnetic systems is quite illuminating. When the coars- 
ening dynamics proceeds, larger and larger magnetized 
domains are present, and the dynamics indeed becomes 
slower and slower. The extent to which slow dynamics is 
always accompanied by a large coherence length in any 
glassy material is an important open problem in this field. 
In this paper, we adopt the poin t of view that this is in- 
deed the case. See Refs. [22], |23( for related discussions. 

Using this spatial description of slow dynamics, 
preasymptotic behaviours were conjectured in Refs. 0,0 
to obey a generalization of Eq. Q. The physical pic- 
ture is the following. At a given wait time, t 2 , of the 
aging dynamics, the system has become locally equili- 
brated up to a coherence length scale, £(t 2 ). It can thus 
be seen as a heterogeneous mosaic of independent and 
quasi-equilibrated sub-systems of size i(t 2 )- Now, a dy- 
namic measurement performed on this mosaic is in fact 
probing an ensemble average over equilibrated systems of 
finite size £(t 2 ). Translating this idea into an equation, 
one arrives at the conjecture that the following general- 
ization of Eq. (@J can be valid 0, 0| : 



X(C,t 2 ) = 



f dq'P{q',L), L = £(t 2 ), 
Jo 



(14) 



which implies that non-equilibrium properties of the sys- 
tem at time t 2 , as encoded in the FDR, are related to the 
equilibrium properties of a system of finite size L = £(t 2 ), 
as encoded in the Parisi function. Note also that static 
and dynamic preasymptotic effects are explicitly present 
in Eq. 1)14(1 . since both finite times, t 2 , and sizes, L, are 
taken into account. 

The conjecture (|14f) was first formulated and tested 
in Ref. |4| in the context of the non-equilibrium critical 
dynamics of the 2d XY model, and further studied in 
Refs. H l24l | in the numerical simulation of the Ising spin 
glass in dimension d = 2. In both cases, the observed be- 
haviours were by construction preasymptotic, since for 
both models the true asymptotic behaviour was simple 
equilibrium with X(C) — 1 and P(q) = S(q). However, 
on numerically accessible time scales, a non-trivial be- 
haviour reminiscent of the one observed, say, in realistic 
spin glasses was found, and the conjectured relation i(Ti|) 
was indeed obeyed to a very good precision. 

It is straightforward to obtain the relation correspond- 
ing to Eq. iTHI) in the case of a driven glassy material: 



X(C,e) 



dq'P(q',L), L = l(e) 



(15) 



where the relevant control parameter in the dynamics is 
now the amplitude of the driving force, e. The advan- 
tage of this formulation, as exemplified in the numerical 
simulations below, is that the dynamics is stationary, so 
that the coherence length itself does not change during 
the dynamical measurement, as is the case in the aging 
regime. 



B. Single-site quantities 

A different generalization of Eq. was recently dis- 
cussed in the literature [JJ, llCj - Inspired by numeri- 
cal results obtained in a disordered spin system 0, the 
argument reviewed in section II Bl was reformulated us- 
ing single-site quantities defined for a given sample |1C| . 
Therefore, disordered systems only are concerned by this 
approach. 

Consider the alternative single-site definition of a gen- 
eralized susceptibility fiof . 



M r 
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E 



ds° 



dh ai . 



(16) 



h=0 



where M thermodynamically coupled copies of the same 
system have been introduced, so that s" denotes the value 
of spin i in the copy a. The role of a small coupling of am- 
plitude k between the M copies is discussed in Ref. [Toj . 
Its presence is necessary to properly define single-site 
overlap distributions Pi(q), although its precise form is 
inessential. 

If an inversion of the limits M — > 00, t — > 00, fc — > 
and /i — > is again allowed, then a reasoning analogous 
to the one developed in section IT*Bl leads to the equality 



Xi{C) 



dq'Pi(q'), 



(17) 



between the single-site FDR and local overlap distribu- 
tions. The single-site FDR is obtained by averaging local 
correlation, C'i(ti,t 2 ) = Si(ti)si(t 2 ), and response func- 
tions, Ri(ti,t 2 ) — dsi(ti) / dhi(t 2 ), over various realiza- 
tions of the thermal history for a single realization of the 
disorder. We note finally that the derivation of this lo- 
cal FDT makes use of the same type of (uncontrolled) 
hypothesis of the global one. 



C. Single-box quantities 

Although Eqs. (|17|) and (|14l) arc different, they are also 
consistent, because Eq. I|17f) physically follows from the 
concept of a local equilibration in space, so that dynam- 
ical properties at site i can be linked to static properties 
at the same site i. 

The disadvantage of this single-site formulation is how- 
ever evident, since it entirely relies on the presence of 
quenched disorder in the system. A slight generalization 
of Eq. I|17fl would be to consider single-box quantities 
instead of single site to get similarly: 



XV(C) 



(18) 



where the box FDR Xf(C) is defined from box dy- 
namical functions, C^{ti,t 2 ) = v^ 1 ^^ Sj(t\)sj{t 2 ) and 
Xi(*i)*2) = v ^ ly E,j 9sj{ti)/dhj{t 2 ), and the sums run 
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over a box of finite volume v centered around site 
i 0, 0> 13- The corresponding static box overlap dis- 
tribution is denned as (25| : 




using the same notations as in J5j. The interest of the 
formulation (|18|l is that it should become independent of 
the considered site i for a sufficiently large box volume 
v. A second interest is that it suggests an alternative to 
Eqs. and l(TS)l . namely 

X(C,t 2 ) = [ dg'P"(g), v=e d (t 2 ), (20) 
Jo 

for the aging case, and 

X(C,e)= [ dq'P v (q'), v = i d {e), (21) 
Jo 

for the driven case. 



III. SIMULATIONS OF AGING AND DRIVEN 
DYNAMICS 

We now turn to a numerical investigation of the two 
relations lTH|) and (|T5j) discussed in the preceding section. 
We first present the model studied and some technical 
details, before giving our results both for a driven and 
an aging dynamics. 

A. Model and details of the simulations 



This implies that conclusions drawn for systems as small 
as L = 10 have to be taken with much caution |27j . 

The aging dynamics is simulated by preparing the sys- 
tem at initial time in a random configuration, thus mim- 
icking an infinite temperature state. The temperature is 
then changed at time t = to its final value. We will 
present below data for the temperature T = 0.7T C in 
both d = 3 and d = 4, where T c = 0.95 and T c = 1.8, re- 
spectively. Dynamical measurements are then performed 
during the resulting aging process. 

To simulate the driven d yna mics, we use the same type 
of perturbation as in Ref. |15|. On each link of the 
lattice, a coupling is added, where Jij is drawn from a 
Gaussian distribution of mean and variance e, which de- 
fines then the amplitude of the driving mechanism. The 
key information about the J's is that they are chosen to 
be antisymmetric, Jy = ~~Jji, which implies that their 
effect cannot be incorporated in a Hamiltonian force of 
the type l|22|l . Hence, the effect of the J's is that of a 
non-conservative force, and the system is thus externally 
driven. As a result, a driven stationary state is reached at 
sufficiently large times, even for temperatures below T c : 
aging is stopped by the driving force [l5|. In that case, 
measurements are performed in the stationary regime, 
meaning that a considerable amount of the data corre- 
sponding to transient behaviours has to be discarded. 

In both cases, dynamical measurements consist of the 
measurements of the global functions and J2J. Since 
these measurements are by now classic, we do not discuss 
them further here [l3|. 

In order to perform a comparison with static data, 
we shall make use of the published equilibrium data of 
Ref. [H in d = 3 and Ref. in d = 4. Fortunately, 
box-overlap data exist in d = 3 and we shall also make 
use of those |25|. 



In this section, we study numerically the Edwards- 
Anderson model of a spin glass defined by the Hamil- 
tonian [2(| 

N 

H = -^J ij s i s h (22) 

where Sj (i = 1 • • • N) are N = L d Ising spins located at 
the sites of a cubic lattice when d = 3, or a hypercu- 
bic lattice when d — 4, of linear size L. The coupling 
constants are random Gaussian variables of mean 
and variance 1 and the sum in Q22[l runs over nearest 
neighbours. We use a standard Monte Carlo algorithm 
where the spins are randomly updated. Times are given 
in Monte Carlo steps, where one step represents N at- 
tempts to update a spin. Sizes are given in units of 
the lattice spacing. To study the dynamics of the sys- 
tem, a large system size, L = 40 (d = 3) and L = 24 
(d = 4), is chosen, so that few realizations of the disorder 
are needed, typically 10. Such large sizes are necessary, 
as the Fig. ^ discussed below will clearly demonstrate. 



B. Time is length 

To accurately test the relations (I14f) and 115fl . the fol- 
lowing steps must be followed. 

• Identify and measure in the dynamics the coherence 
length i. 

• Measure global correlation and response functions, 
from which the corresponding FD plot is built. 

• Measure in equilibrium simulations the overlap dis- 
tribution functions. 

• Compare static and equilibrium data with carefully 
chosen values of the parameters (size, time, driving 
force), as suggested by Eqs. {13} and (JTSJ. 

The first point, the identification of a coherence length, 
has been discussed for the aging dynamics of the Ising 
spin glass in various dimensions in Refs. [3(], On], O^] , so 
that the existence, temperature and wait time behaviour 
of this quantity is well-characterized. 
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FIG. 1: The correlation function $23} for d = 3, T = 0.7T C 
and various amplitudes e of the driving force. 

To the best of our knowledge, no such results are avail- 
able for the driven dynamics. We have thus measured in 
the driven stationary state the standard four-spin corre- 
lation function 

i N 

C 4« = ^E< s i<+r s i s U' (23) 
i=i 

where the overline means an average over the disorder, 
and (a, b) are two independent copies of the system. 
Some corresponding curves are shown in Fig.Q for d = 3. 
Very similar results are also obtained for d = 4. The fea- 
tures observed in this figure are expected. The amplitude 
of the driving force is known to control the relaxation 
time, tj-p] = t re i(e), of the system, the smaller e the larger 
^rei |l5|. From Fig.^ we recognize that a slower dynam- 
ics also implies the existence of a larger coherence length, 
£(e), as revealed by a slower spatial decay of C±(r). 

A non-ambiguous definition of the coherence length is 
obtained, as in the aging regime, through the study of 
the scaling properties of Ci(r). For the aging, it is known 
that the correlation l|23[l is well described by the scaling 

form mm 

which can be seen as the definition of the coherence 
length I. We find that the same scaling behaviour 124|) 
is obeyed in the driven dynamics as well, with the same 
value of the temperature-dependent exponent a(T) as in 
the aging regime. As an example, we report the scaled 
data for T = 0.7T C and d = 3 in Fig. H We find that 
the scaling function C(x) is well described by a simple 
exponential form, C[x) = exp(— x). This contrasts with 
the 'compressed' exponential form reported in the aging 
regime [33. We have no explanation for this difference, 
which is however inessential. 

The results of Figs. and [5] allow us to obtain the 
relationship between the coherence length t and the am- 
plitude of the driving force, e. From the time depen- 
dence of the correlator it is also possible to extract 



FIG. 2: The data of Fig.^Jrescaled according to Eq. pijl with 
a = 0.5, T = 0.7T C , d = 3. The straight line indicates a simple 
exponential behaviour of the scaling function C(x) — exp(— x). 

a relaxation time, i ro i(e). It is thus natural to eliminate 
e to obtain the relationship between time and length. 
Not surprisingly, we find that the relation £(i re i) in the 
driven regime, is well compatible with the relation Ifa) 
measured in the aging regime at the same temperature. 

Although this coincidence might seem anecdotic in this 
context, we emphasize its physical importance. This sug- 
gests indeed that whatever the control parameter for the 
dynamics is, waiting time or driving force, the dynam- 
ics is such that the same dynamic scaling £(t) is obeyed: 
'Time is length'. 

C. Generalized FDT for the driven dynamics 



Following the program presented in section IIII Bl we 
compute in the driven dynamics correlations and suscep- 
tibilities, and build the corresponding FD plots. Collect- 
ing the static data of Refs. [23. |2§1 l29j , we are in a posi- 
tion to start and investigate the validity of the relations 
d, (E>|, J!]) and HI). 

We present in Fig. |3 the results obtained in d = 4. In 
this figure, there are three different FD plots, represented 
with points, and corresponding to three different values 
of the driving force, and thus to three different coherence 
lengths £(e). The lines correspond to the double integral 
S(C, L) defined in Eq. O of the Parisi function P(q, L) 
for three different system sizes. The agreement between 
the two sets of curves is extremely good. That this is 
not a mere coincidence is supported by three important 
facts. 

• The two types of measurements are performed en- 
tirely independently, in completely different physi- 
cal situations; 

• The agreement between the two sets of curves is 
very good on the whole range of the curve, not on 
a single point. If one remembers that each curve is 
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FIG. 3: Test of the static - dynamic relationship (1151 for the 
driven dynamics for d — 4 and T — 0.7T C . The points are 
dynamic FD plots, the lines are the double integral of the 
Parisi function measured in equilibrium simulations. 



built from a rather complex manipulation of vari- 
ous quantities without any fitting procedure or free 
parameter, the agreement becomes much more im- 
pressive; 

• The coherence lengths measured in the dynamics 
are £(0.2) = 1.55, £(0.16) = 2.0, £(0.12) = 2.5. 
They are thus in ratio 1:1.29:1.61. The correspond- 
ing static length scales are in ratio 1:1.33:1.66. 
The agreement with the prediction (|15fl that both 
lengths should scale similarly is thus clearly excel- 
lent. 



D. Another observable: The link overlap 

To test further the hypothesis of a strong link between 
off-equilibrium properties and static ones resulting from 
the local equilibration of the material, we focus in this 
section on a second pair of dynamical observables. Phys- 
ically, this is motivated by the picture given in section 
III Al where an off-equilibrium system was described as a 
mosaic of equilibrated sub-systems of size given by the 
coherence length. A logical consequence is that the link 
between static and dynamic quantities should exist for 
all physical observables. This is equivalent to the predic- 
tion and test of the unicity of the effective temperature 
in structural glasses 0> El EE ■ 

To check this hypothesis, we have investigated the fol- 
lowing quantities: 



N 2d 



2Nd 



N 2d ' 



= ^EE 



dsi(ti)sj(ti) 



(26) 



FIG. 4: Test of the static - dynamic relationship list for 
the driven dynamics for d = 4 and T — 0.7T C for the link 
observables. The points are dynamic FD plots, the lines are 
the double integral of the link overlap distribution function 
measured in equilibrium simulations. Note the coincidence 
between the values of e and L with the ones reported in Fig. 01 



where the prime in the second sum indicates that it is 
restricted to the nearest neighbours of the site i. The field 
hij is thermodynamically conjugated to the observable 
SiSj, and can thus be seen as a random perturbation 
of the coupling constants. From these two dynamical 
functions, a new FD plot can be built. 

Correspondingly, the static quantity to use is the link- 
overlap distribution function, P(qi), defined by 



N 2d 



lim ( 6 , 

Af^oo \ \ 2Nd 



— W s a s a s l3 s f: 



(11 



1=1 j=l 



h=0 



(27) 

We use here the static data published in Ref. 29] 

The comparison between the new sets of static and 
dynamic quantities is reported in Fig. for the same 
parameters as in Fig. El It should not come as a surprise 
that the range of correlators covered by Fig. 0] is smaller 
than the one covered in Fig. The long-time limit of 
the correlator is indeed non-zero in that case. 

It is again clear that the agreement between static and 
dynamic data is very good. The point to emphasize is 
that for a given amplitude of the driving force, the same 
size L as in Fig. [3] is used to compute static data, in 
agreement with the physical picture described above. 

These results are a striking confirmation of the local 
equilibration hypothesis formulated in this paper. 



E. Using the box overlap distributions 

Having established numerically the validity of the rela- 
tionship between static and dynamic properties in d = 4, 
we now turn to d = 3. This will allow us to make 
use of the box-overlap distribution functions published 
in Ref. pEj . To compute this function, the linear size of 
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FIG. 5: Test of the static - dynamic relationship 12111 for 
the driven dynamics for d = 3 and T — 0.7T C . The points 
are dynamic FD plots, the lines are the double integral of 
the box-overlap distribution function function measured in 
equilibrium simulations, B refers to the linear size of the box, 
L to the linear size of the total system. 



FIG. 6: Test of the static - dynamic relationship 1201 for 
the aging dynamics, for d = 3 and T — 0.7T C . The points 
are dynamic FD plots, the lines are the double integral of 
the box overlap distributions measured at equilibrium. The 
discrepancies for the short waiting time ti = 27 are discussed 
in the text. 



the total system was fixed to L = 12, and the linear size 
of the boxes, B, was varied [25| . 

Our numerical results are reproduced in Fig. El where 
the FD plots built from the spin-spin dynamical functions 
JTJ and © for various amplitudes of the driving force 
are compared to the static data obtained from the box- 
overlap distributions. 

As for d = 4, we find a very good agreement between 
static and dynamic data. The coherence length scales are 
in ratio 1:1.35:1.59:1.88, while the size of the boxes are in 
ratio 1:1.33:1.66:2, the coincidence being again satisfying. 

We conclude that numerical simulations in d = 3 and 
d = 4 of the driven dynamics of the Ising spin glass nicely 
support the theoretical propositions (|15|) and l|21|) . 



F. Generalized FDT in aging dynamics 

We conclude the report of our numerical results with 
a brief mention of the results we have obtained in the 
aging regime for d = 3,4. Similar results have already 
been published, as mentioned above, but the crucial co- 
incidence between time (in the dynamics) and size (in the 
statics) was not noted, although it was supported by the 
data [13, • We have also discussed above the results 
obtained for d = 2, when no spin glass phase exists 0,13- 

We give an example of our own data in the aging regime 
in d = 3 in Fig. El There are four important points we 
want to make on these data. 

First, we make use in this figure of the box-overlap 
distribution function, which was not the case in previ- 
ous works. The total overlap distribution functions give 
however results which are similarly good. 

Second, the dynamic data can be described using the 
same equilibrium quantity P(q, L) as for the driven dy- 



namics. This is a clear support of the analytical result 
that driven and aging dynamics should give rise to similar 
deviations from the FDT p^lll[l5|. 

Third, we have intentionally included data for a very 
short wait time, — 27, for which the agreement with 
Eq. I|20|) is not quite as good, to underline the fact that 
during the dynamical measurements the coherence length 
itself evolves with time. This was obviously not the case 
when the driven dynamics was considered. For this short 
wait time, for instance, the final time is t\ = 10 5 , so that 
we do not expect the FD plot to follow the line built from 
static data obtained at a given box size corresponding to 
the earlier time L = ifo)- Indeed, a clear deviation at 
large times (corresponding to low values of the correlator) 
is observed. This was anticipated and noted in Refs. 0, 
5] , and was another motivation to focus also on the driven 
dynamics in this paper. 

Fourth, we have discussed above the link between the 
coherence length £ and the system size or box size in 
terms of relative ratio. It is clear that an absolute com- 
parison is not possible, since the coherence length has 
no 'absolute value', contrary to the box-size univocally 
defined in the static computation. However, the 'large' 
numerical factor between L and I reported in Ref. |2Cj | 
(about 4 in d = 3) was an argument used against the 
theoretical proposition (|14|) formulated in this and previ- 
ous works [33. It can be seen from Fig.^that the factor 
between I and B is less than 2 when the box overlap dis- 
tribution is used, thus clearly weakening this criticism. 

IV. CONCLUSION 

In this paper, we have given a scaling argument to sup- 
port the proposition of Refs. to generalize the rela- 
tionship between the off-equilibrium FDR and the equi- 
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librium overlap distribution to the physically relevant do- 
main of finite times and sizes. Inspired by previous ap- 
proaches to the asymptotic relation [2, H3 > we proposed 
the relations i|14|) and (|15fl between the FDR computed 
at any moment of the dynamics, and the equilibrium box 
overlap distribution function, where the size of the box 
to be used is imposed by the dynamical coherence length. 

We have argued that these relations reflect the sim- 
pler notions of the local equilibration in space and the 
dynamical heterogeneity of the system, two ingredients 
which arc inherent to systems with slow dynamics, inde- 
pendently of their space dimensionality, their dynamics, 
or the presence of disorder. 

We have emphasized throughout the paper that the 
relations l|14|l - <|15|) were the physically relevant ones, as 
compared to the asymptotic form of Eq. (£Q). Hence, they 
have to be taken into account before drawing conclusions 
on the equilibrium glass phase one seeks to study |20j| . 
This remark, already made in Refs. 0, 0, is repeated 
here because there recently appeared FD plots built from 
experimental dynamic measurements in spin glasses |l8j . 
As anticipated Q , this dynamical measurements are very 
similar to the ones observed in simulations, despite the 
enormous difference in time scales between experiments 
and simulations. This is because the growth of the coher- 
ence length is so slow in spin glasses that this difference 
in time scales reduces to a very modest factor in terms 
of the coherence length pjj . I32L |3S| . This in fact implies 
that the 'finite size effects' reported in simulations, are 
present in experiments as well. This points towards the 
experimental irrelevance of the thermodynamic limit in 
spin glasses. 



Finally, we made clear that the relations l|14|) - l|15(l 

had not the character of a theorem, meaning that more 
work is needed to assess their range of validity. In par- 
ticular, coarsening systems are a simple counterexample 
of this relation. This is due to the fact that the dynami- 
cal behaviour is entirely dominated by topological defects 
which are absent in the equilibrium simulations, so that 
finite time dynamics and statics do not coincide |3.l5ll39|. 
There remains to be seen to what extent this alternative 
'defect' description of glassy dynamics generally holds in 
systems like supercooled liquids or soft glassy materials. 
This is certainly a quite challenging issue. 
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